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Abstract 

The supersymmetric approach is used to analyse a class of two-dimensional quantum 
systems with periodic potentials. In particular, the method of SUSY-separation of 
variables allowed us to find a part of the energy spectra and the corresponding wave 
functions (partial solvability) for several models. These models are not amenable to 
conventional separation of variables, and they can be considered as two-dimensional 
generalizations of Lame, associated Lame, and trigonometric Razavy potentials. All 
these models have the symmetry operators of fourth order in momenta, and one of 
them (the Lame potential) obeys the property of self-isospectrality. 
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1 Introduction 

The supersymmetric (SUSY) approach has provided a powerful impulse for 
new developments in analytical studies in Quantum Mechanics. To date, most 
new results have been obtained for one-dimensional quantum systems (see the 
monographs and reviews [1]). Over the last two decades, indubitable progress 
has been achieved also for higher-dimensional systems within the framework 
of SUSY Quantum Mechanics [2], [3], [1], [S]. The most essential results have 
been found for two-dimensional Quantum Mechanics. In particular, the new 
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approach for the construction of completely integrable systems with symmetry 
operators of fourth order in momenta was proposed [3]. Two new methods - 
SUSY-separation of variables and two-dimensional shape- invariance [S], |1] 
- have allowed the problem of solvability to be tackled for two-dimensional 
quantum systems, beyond the standard separation of variables. 

The method of SUSY-separation of variables has been applied successfully to 
investigate the spectra and wave functions of some models on the whole plane 
X = {xi,X2), which are not amenable to conventional separation of variables: 
the Morse potential [5], [6], the Poschl- Teller potential [7], [8], and some 
others [H], i). 

The starting point of this approach is the analysis of solutions of the in- 
tertwining relations between a pair of two-dimensional Hamiltonians of the 
Schrodinger form: 

Hix)Q- = Q-H{x)- Q+H{x) = H{x)Q+- (1) 
i7 = -A(2) + V(x); = -A(2) + V(f). (2) 



In general, the intertwining operators are the operators of second order 
in momenta, i.e. in derivatives. In [7j (see also [H]), the particular class of 
intertwining operators was considered: operators with twisted reducibility 



Q- = {Q 



-di + dixix)){a3)iki+dk + dkxix)), di 



d_ 
dxi 



(3) 



where x{^)iX{^) are two different functions (superpotentials), is the Pauli 
matrix, and summation over /c, / = 1, 2 is implied. In this class of models, both 
Hamiltonians H and H are quasifactorized [7] 



H = i-di + dixim+di + dam = -A(') + mm' - mxix), 

and hence their energy spectra are non-negative. 



(4) 



It was shown in [7j, [8] in a general form that Eqs.Q and ([S]) lead to the 
following representation for superpotentials X; X terms of four functions 
/^i,2,± : 

X = /ii(xi) + iJ,2{x2) + /i+(a;+) + /i_(a;_), 
X = fJ-iixi) + iJ,2ix2) - - 
where x± = {xi ± X2)/V^- These functions /ii,2,± satisfy the equation 



fJ^'liXl) + IJ,2{X2) - 



0. 
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By = /i' substitutions, this becomes a purely functional equation with no 
derivatives: 

[(f)+{x+) + (f)-{x-)] = -Mx2) - 0-(a;-)] • (5) 

The solution of Eq.(|5]) is necessary to build the solutions of intertwining rela- 
tions for the potentials V, V and for the supercharges Q^- In particular, 

V{x) = (0?(xi) - <f>[{x^)^ + (^4{X2) - <l>',{x,)^ + (0^(x+) - 0'+(a;+)) + 

+ (02_(x_)-0'_(x_)), (6) 
Q± =0^ - dl ± v^(^0+(x+) + T V2(0+(x+) - 0-(x_))92 - 

+ {<j)\{x2) -02(a;2)) +20+(x+)0_(x_). 

Thus, in order to find the systems with intertwining ([T]) by supercharges of the 
form ([3]), it is necessary to solve ([s]). This equation seems to be rather compli- 
cated, but it appeared to be solvable in a general form. The two-dimensional 
generalization of the Poschl- Teller potential, investigated in [7j, [8], was based 
just on a particular solution of Eq.([5]). 

In the present paper we focus our attention on the 0i,2,±; solutions which 
are periodic in the variables X\ and X2- These solutions will be further used to 
build a class of two-dimensional potentials V{x) of the form ([6]) - not amenable 
to standard separation of variables - which are periodic along X\^X2 with the 
same periods. By means of the method of supersymmetric separation of vari- 
ables [3], |1] we shall derive the partial solvability of these periodic models: 
several energy eigenvalues and corresponding eigenfunctions will be found. 
Until now the supersymmetric approach has been used for the analysis of pe- 
riodic potentials in one-dimensional case only (see for example [lOj, [IT], 
[12]). To the best of our knowledge, this i^the first attempt to study ana- 
lytically two-dimensional periodic potentials, not amenable to separation 
of variables, within the framework of SUSY Quantum Mechanics. 

It is appropriate to make some remarks concerning the possible spectra of two- 
dimensional periodic systems. The general statements are well known from 
textbooks (mainly on solid state theory) [14j: the spectra of rf— dimensional 
[d > 2) models with periodic potentials in general have a band structure sim- 
ilar to that of the one-dimensional case. However, in contrast to d = 1, no 



See some remarks in [13] (Section 5) and [8]. 
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strict results on the (anti)periodicity properties of the band edge wave func- 
tions are known [T^, and therefore analysis of the spectra of two-dimensional 
systems is much more complicated. In some sense, the situation resembles 
the non-periodic case: no analogues of the oscillation theorem are known for 
d > 2 quantum systems. This is a reason why the analysis of multidimen- 
sional excited bound states is much more difficult than in the one-dimensional 
situation. 



To imagine the variety of possible structures of the band spectrum, one can 
consider (contrary to the rest of this paper) the simplest two-dimensional 
periodic systems which allow separation of variables (see, for example, [TB] 
and references therein). After separation of variables, both one-dimensional 
problems have the band structure of energies ei(fci), e2{k2). Let us assume that 
they have a finite number (one or two) of band gaps (as in Subsections 4.2 
and 4.3 below). Then, the positions of the two-dimensional band edges depend 
crucially on the parameters of one-dimensional bands ei, e2. In particular, these 
bands may be overlapped, so that the band gaps (or at least some part of 
them) of the two-dimensional spectrum, E = ei + 62, may even disappear. It is 
natural to consider two limiting cases: of almost free particles (the band gaps 
are vanishing) and of tight binding particles ( the band gaps are very wide). 



The structure of this paper is as follows. Section 2 is devoted to the solution 
of the functional equation ([s]): the explicit expressions for all four functions 
0i,2,± will be found. In Section 3 the symmetry properties of the general so- 
lution of Eq.([5]) are analysed. In Section 4, which is the main part of the 
paper, a new class of two-dimensional systems, partially solvable by means 
of SUSY-separation of variables, is constructed. The potentials of these sys- 
tems are periodic on the plane, and they can be considered as two-dimensional 
generalizations (not amenable to separation of variables) of fairly well known 
[T7] . [To] . pT] . [12], [IB], [I9] one-dimensional potentials: the Lame potential, 
the associated Lame potential and the trigonometric Razavy potential. Using 
the known positions of band edges and corresponding wave functions for these 
one-dimensional models with a finite number of gaps, some energy eigenvalues 
and their eigenfunctions for two-dimensional periodic generalizations will be 
obtained exphcitly (partial solvabihty). Section 5 includes a discussion of cer- 
tain specific properties of the models of the previous Section and some limiting 
cases. In particular, all two-dimensional models constructed are integrable: the 
symmetry operators of fourth order in momenta commute with the Hamilto- 
nians. In addition, the two-dimensional Lame model obeys the property of 
self-isospectrality, like its one- dimensional prototype [TU] . 
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2 The general solution of the functional equation 



In this Section we analyze the functional equation (jsj), which plays an im- 
portant role in our approach. This analysis was started in the Appendix of 
the paper [7], where the necessary conditions for the existence of its solutions 
were derived. First, we will remind briefly the main steps of this derivation. In 
the formulas below we shall imply that the functions (j)i^2, 4>± depend on the 
corresponding arguments: 0i,2 = </'i,2 (2^1,2) and 0± = 4>±{x±), unless otherwise 
stated. 

Acting by {Of — on both sides of (j5|, we obtain: 

2(^(^1^ 0l9i-^0'252)(0--0+)= 01^(0^-0+). (7) 

It is now convenient to introduce a new unknown function, A : 

0_ - 0+ = 01 101021"^^^ A(a;i, X2). 
Substitution of this definition into ([T]) gives: 

( r^. + r*) A = 0. 

\0i 02 / 

and therefore A depends only on J^^ ^/^^^ ~ W^^^^ ' ^^"^ general 
solution of ([7]) is: 

From the initial equation ([s]) we also have that: 



From Eqs. ([s]), (|9| we could already obtain the solutions for 0-|-, but we have 
to check that these solutions will indeed depend on proper arguments. Thus, 
the constraints for the function A and functions 0i_2 are obtained from the 
equations: 

9±0^ = 0. 

The result is: 
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(10) 

(11) 



Here we disregard the trivial solution A = 0, for which 0+ = 0_ 
arbitrary, but the potentials ^ are amenable to separation of variables. 



0, 01,2 are 



Otherwise one can exclude A, dividing Eq.(lO) by Eq.(ll): 



202 - 201 + 0101 - 0202. 



(12) 



There is no separation of variables in (12), but it will appear after applying 
the operator 8182.1 so that: 



,2^/ 



2a = const. 



Integrating, multiplying by 0^ 2, integrating again and taking into account 
(12), one has that 0i 2 must satisfy the equation: 



k'\2 



1,2 = «01,2 + ^01,2 + C, 



(13) 



where a, b, c are arbitrary real constants. All solutions of this equation can be 
expressed in terms of elliptic functions, and they are described for different 
ranges of parameters, for example, in [20] . 



The subsequent discussion will depend crucially on the sign of the discriminant 
D = ~ 4ac. In the present paper we restrict ourselves to the case of D > 0, 
i.e. when the quadratic polynomial 00^ 2 + ^0i 2 + ^ has two different real roots, 
which will be denoted as ri and r2 (let ri > r2). For positive values of the 



discriminant, three types of solutions of (13) exist, depending on the relative 



position of the roots ri,2. These solutions are proportional either to '^"■^"^l^) 



tod^n(ax\k) . 1 
7 'ryi or to — f i,s 

sn(ax\k) sn(ax\k) 

functions [TT] . 



sn{ax\k) 

, where sn,cn, dn are the well-known elliptic Jacobi 



Equations ( |I3| ) are the necessary conditions for the functions 0i,2 to satisfy 
equation (|5|). But are these conditions also sufficient? To answer this question. 



■^1,2, 



we must solve Eqs.(lO), (11) for function A, with arbitrary elliptic functions 



satisfying (13). From Eq.(13), the argument of the function A can be 



written in the form: 



M^,. _ 02 

J 0UO ^ J 0'2 



iv) 



dr] 



1 



2a(ri - r2 



-In 



r2 



(14) 
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The module sign is not important here, since 



u = 



> 



(15) 



because of (13). In terms of this new variable Eq.(lO) takes the form 

K^) ~ ~4z/(z/i - a)' 
where the notation a = sign{(f)[{xi)(t>2{x2)) was used. 
The solution of this differential equation is 



Z/4 



A(i^) = Ao — 

1/2 — cr 



Ao = const. 



Formally substituting (17) into the l.h.s. of Eq.(16), we obtain: 



z/2 + cr — Avd^a 
4z/(z/^ — cr) 



(16) 



(17) 



Let us prove that the extra term in the nominator of the r.h.s. actually van- 
ishes; i.e. that vd^a = 0. The variable u may be considered as a function of 
(pl and and hence 



d d(j)\ d 



dv dv d(f)i 



+ 



hi d 



dv d(j)2' 



In turn, due to Eq.(|5]), one can rewrite derivatives over 0^ in terms of deriva- 
tives over d)'- : 



d 



A2 PiA' 



a{2(t)j - ri - rz) 



--sign[ 



2^a(02-rO(02-r2) 
Gathering everything together and using = 25(0'^), we find that: 

which is zero for all values of Xi^2 since at the points where 5(0j) ^ 0, the 
denominator has no singularity (it tends to ±(ri — r2) for (p'- 0). 



^ This solution coincides with the solution in the simpler case of a constant (not 
depending on z/) value of a. 
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Now, according to Eqs.(|8|, ^ solution (17) gives 0± in terms of functions 

01,2: 

'>2 T 01 



Ao- 



1^ 



2 Z/2 — (7 



The important fact is that without any new constraints on 0i,2(a;i,2) (besides 



Eq.(13)) the functions 0± depend on the proper arguments x± : 



For any pair of functions 0i,2 satisfying (13) there exists a pair of functions 
0± given by (IS), such that they are solutions of 



Closing this Section we shall derive an alternative expression for (p^, to be 
used later (see Section 4). Let us define the function 



G 



+ 



+ h(t)i + h(t)i + 2c 



One can check by straightforward calculations, that 



Then, taking into account 



(19) 



Eq.(19) can be rewritten in compact form as: 

d 



G 



0. 



Hence, 0?. = Const ■ G, where the constant can be calculated by considering 
this expression in some specific point: Const = Ag/ (4(6^— 4ac)), with arbitrary 
constant Aq. We therefore obtain the desired alternative expression for (f)_ : 



Al 20;0'2 + 2a(pl(pl + 60? + b(f)l + 2c 



b'^ — Aac 



(20) 



By a completely analogous calculation, or directly from Eq.(|5]), the expression 
for 0+ can be derived: 



A^ 



6^ — 4ac 



20'i0'2 + 2a0f0^ + b(f)( + b(j)^ + 2c 
(02 + 0i)2 ■ 



(21) 



Formulas (20) and (21) are very convenient to prove that 0± satisfy the same 



8 



Eq.(13), but with different coefficients. Direct calculations show that: 



4: 



2 

^^0 



2 {b^ - 4ac)4 - 



^1 + ^ 



(22) 



Let us note that although the discriminant D = IP' — 4ac was chosen to be 
positive, the analogous discriminant D for Eq.(22), which is equal to D = Aac, 
can be non-positive as well. 



3 Symmetries of the functional equation 



As pointed out in a previous paper [7] (Subsection 3.2), the functional equation 
(|5| has two discrete symmetries: 5*1, 5*2. After the change of variables yi^2 = x±, 
y± = [yi ± y-2)l\pl = Xi,2, equation becomes: 



My+M+ivi) + 0-(y2)] = -My-M+ivi) - </'-(i/2)], 

which, by rearrangement of terms, can be brought to: 

+ My~)] = -(p-{y2)[My+) - My-)]- 

The last form coincides with ([5| up to interchanged (^1^2 and Hence, if 
(01, 02, 0+, 0_) is a solution of |5]), then the set (0+, 0_, 0i, ^2) is also a solu- 
tion (this discrete symmetry was called 5*1). By writing the four functions in 
brackets, we mean that the first of them should be put in the place of 0i in 
(|5|, the second one in the place of 02, and the last two in the place of 0+ and 
0_, correspondingly. 



This symmetry can be observed explicitly from the general solutions (20) 



(21). Namely, if we calculate 0i,2 from (f)± with the use of analogues of (20) 



(21): 



~2 Ag 20V0-_ + 250^0^. + b<Pl + b<Pl + 2c 

"^''''b^-A-ac' (0-±0+)2 ^ ' 

then we must arrive at the same functions 0i,2- Indeed, calculating the r.h.s. 



~ 2 ~ 

of (|23[), one obtains 0i = Ao0i/(8c), and by proper choice of the arbitrary 



constant Aq (which can be imaginary) one has 0i = 0i. Thus, 02(2/2) = 
0i(?/i)(0+(?/+J + 0_(y_))/(0_(i/_) - 0+ (?/+)), and by comparing it with (g 
we see that 02 = 02 too. 

There is another symmetry S'2. If (0i, 02, 0+, 0-) is a solution of ([S]), then 
(01, —02, ^) is also a solution of (5). This was shown in ^ directly from 
the functional equation ([s]), and now we shall derive it from the general solu- 
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tion (18 



). It can be rewritten in the form: 



Ao 



z/2 + cr 



1 

z/2 



O" 



If one performs the change 



An 



^1, 



92 =F 01 

bo ± 



then (7 



-0", and 



y\b\ 

which completes the proof. 



z/2 — cr 



Z/2 + (T 



-02 T 



-09 ± 



Const 



It should be noted that these symmetries may play an important role. In 
particular, the first one -Si - was applied in |7j to combine the shape-invariance 
with SUSY-separation of variables in the 2D Poschl- Teller model. 



4 SUSY-separation of variables for some models with periodic po- 
tentials 



4-1 The algorithm of SUSY-separation 



In the general expressions ([6]) for the potential V{x) one sees the very special 
dependence on coordinates {xi,X2): there are two terms that do not mix xi 
and X2, and two mixing terms that depend on the "light-cone" variables x± : 

V{X) = Vi(Xi) + V2{X2) + V+{X+) + V1(X_); Vi,2,± = ^l2,± - 0'l,2,±, (24) 

all terms being represented in " supersymmetric" form. Just the last two terms 
of V{x) show that separation of variables in the Schrodinger equation is not 
possible. 

The method of SUSY-separation of variables (see details in [S], [7j) allows 
us to partially solve some two-dimensional models, i.e. to find a part of their 
spectra and corresponding wave functions. Due to the intertwining relations 
Q, the subspace spanned by zero modes f2n(a?) is closed under the action 
of the Hamiltonian H. Therefore, the wave functions \E'(a;) can be built as 
linear combinations of zero modes ^ni^) of the supercharge Q~^. In turn, 
the zero modes flni^) can be found through the specific similarity ("gauge") 
transformation of the supercharge ([s]) to the separated form: 

q+ = e-^^^^Q^e^^'^ = df - 9| - (0? - 0',) + (0^ - 0'^) , (25) 
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where the function k is given by: 



K[X) = — 



(26) 



This transformation allows us to reduce the two-dimensional quantum problem 
for Qn{x) to a pair of one-dimensional problems: 



where pn and rjn are the eigenfunctions of one-dimensional Schrodinger equa- 
tions: 

-Vn + i<Pl - <P'l)Vn = enVn, (27) 
- Pn + (02 - 02)Pn = ^npn (28) 

[en are the constants of separation). Let us note that the one-dimensional 



potentials in (27), (28) coincide exactly with the first two terms in the two- 



dimensional potential (24). Moreover, they take the form typical of potentials 



within the framework of one-dimensional SUSY Quantum Mechanics [T]. In 



this sense two-dimensional potentials (24) can be thought of as generaliza- 



tions of one-dimensional models (27), (28). 



Our goal in this Section is to build a new class of partially solvable models 
within the framework of the approach described above. In practice, this means 
that we have to present a new class of functions 0i,2 (3^1,2) , which obey the 
following special requirements: 



they must obey the functional equation ([s]), i.e. Eq.(13). 

the one-dimensional models (27) and (28) with superpotentials 0i^2 must be 



exactly solvable in order to provide for the normalizable zero modes fin- 



In such a way we hope to find (j)± by means of (18) (or by means of (20), (pi])), 
and hence to build the two-dimensional potential (|6]). The model with this 
potential will be partially solvable by means of SUSY-separation of variables. 



At the last stage, one has to build the "gauge-transformed" Hamiltonian: 



h{x) 



-k(x)^| 



x)e 



k{x) 



dl+ 



(29) 



25 



and, by direct calculations, find the matrix C*, such that 

huo = Cuj; uj = {ujq, uJi, ... ,ujn)- 
This can be done by the action of the operator h on zero modes a;„ : 
hun = [2e„ + 72(0+ + 0_)9i + 72(0+ - <P-)d2]u;n 



(30) 
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after expressing the r.h.s. linear combination of Uk- 



According to the prescriptions of SUSY-separation of variables (see more de- 
tails in [S], [7j), the eigenvalues of the matrix C give part of the spec- 
trum of the Hamiltonian H. The corresponding wave functions '^k{S) are ob- 
tained as: 

#fc(x) = m{x), 

where the matrix i? is a solution of the matrix equation 

BC = Kb (31) 
with the diagonal matrix A. Its diagonal elements coincide with E^. 
Up to now, this program was performed successfully for the following cases. 

1) A two-dimensional generalization of the Morse potential [5], |lj with 

01 (x) = 02 (a;) ~ exp {ax). 

2) A two-dimensional generalization of the Poschl- Teller potential [7j, [8] with 

01 (x) = — 02(a;) ~ sinh^Vax). 



In both models, the functions (f)±, which were originally guessed directly from 
Eq.([5]), can now be obtained by means of Eqs.(18), (20), (21). 



Below we shall consider new specific choices of functions 0i_2 that will satisfy 
all the aforementioned conditions. The common feature of these new models 
is that the functions 0i,2 are periodic and therefore lead to partially solvable 
two-dimensional potentials, V{x), periodic on the plane (a;i,a;2). 



4-2 The two-dimensional Lame potential 



From the variety of elliptic functions 0i2, which give solutions of Eq.(13), 



we have to choose the subclass leading to exactly solvable one-dimensional 



Schrodinger equations (27) and (28) with periodic potentials. Let us start 



from the particular 2K-periodic solution^ 



0l(x) = 02(x) 



2sn{x\k)cn{x\k) 



dn{x\k) 



(32) 



One can check straightforwardly that this expression satisfies Eq.(13). It can be 



transformed to the third type of solutions of (13), mentioned in Section 2, by the 



Landen transformation |17j accompanied by a suitable shift of argument. 
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where the standard notations [T7] for Jacobi eUiptic functions sn{x\k), cn{x\k), 
dn{x\k) with modulus k G (0,1] and the complete elliptic integral K{k) = 
sin 6 were used. Later on, for simplicity we skip the argument 
k as long as no confusion appears. 



Thus, the one-dimensional potentials Vi^2 of (24), which coincide with poten- 



tials in (27) and (28), have the form of the Lame equation: 



^1,2(0;) = 2esn\x) - k^ 



(33) 



We notice that (33) is only the simplest / 
potential V = l{l + l)k'^sn'^{x) 
I > 1, the superpotentials 



1 case of the general Lame 
with arbitrary integer /. For higher values. 



^2 do not satisfy the basic Eq.(13), and for this 



reason are not considered here. Since the elliptic function sn(x) is 4i^-periodic, 
and 2/('-antiperiodic sn{x + 2K) ■ 
the period 2K. 



-sn(x), the potential (33) is periodic with 



The spectrum and the wave functions for the Lame potential (33) are known 
exactly [17J: it has one bound band with energies e G (0, 1 - k^) and the 
continuous band with energies e G (1, 00), with the band gap between them. 
The band edge eigenfunctions are 



eo =0; 
ei =1 - 

€2 =1; 



i^o^x) = dn{x); 
ipii^x) = cn{x); 
il)2{x) = sn{x), 



(34) 



and the wave functions obey the well known pLTj (anti)periodicity property 
of band edge eigenfunctions for one-dimensional periodic potentials: V'o,i(^ + 
2K) = ^0,1 (x), + 2K) = -iJ2{x). 



According to formulas (20), (21), the choice of (32) for 
compact explicit expressions for 0±(a;±) : 



,2 (2^1,2) leads to very 



^+{x) = -0_(x) = B 



cn 



X) 



sn 



X) 



(35) 



where the new constant B = Ao/(4Vl - F). 

The two-dimensional potential V{x) can be obtained from ([6]): 

V{x) =2k^{sn^{xi) + sn^{x2) - 1) + 

^B{B + V2dn{xi + X2)) , B{B - V2dn{xi - X2)) 



2B' 



(36) 



sn'^{xi + X2) sn2(xi — X2) 

This can be considered as a two-dimensional generalization (not amenable 



to separation of variables) of the Lame potential (33). The potential (36) is 
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a) 



b) 




d) 



3.42. 



0' '^2 (^^^ 



Fig. 1. Panel a): The two-dimensional Lame potential (36) for the parameters 
k'^ = 0.30 and B = 0.05. The period of potential along xi,X2 is 2K{k) 
We used a cut-off of the potential at ±2.0. 

Panels b), c), d): The squares of three known wave functions ^'q, ^'^ 
Eqs.(38)) for two-dimensional Lame potential with the same values of parameters 
and suitable cut-offs. The energy eigenvalues are Eq = 0, Ei = 1.38, E2 = 2.02. 

periodic under the shifts of xi^2 with the periods 2K. The coefficients of both 
terms in (36), singular at (xi ± X2) 0, are such that no fall to the center 
occurs for arbitrary value of parameter B. The plot of the two-dimensional 



potential (36) is represented in Fig. la for a specific choice of parameter values. 



To start the procedure of SUSY-separation of variables, it is easy to calculate 



the function k{x) by (26): 



B 
2V2 



In 



[1 - dn{xi - X2))(l + dn{xi + a;2))' 
[l + dn{xi - X2)){1 - dn{xi + X2)) ^ 



(37) 



The functions Un for the one-dimensional energies e„ at band edges are: 
uJo{x) = dn{xi)dn{x2); uji{x) = cn{xi)cn{x2); uj2{x) = sn{xi)sn{x2)- 



We must act through the operator h (see (29), (30)) on all ui^ in order to find 
the matrix C and its eigenvalues. For this model, the action of the "gauge- 
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transformed" Hamiltonian h on Ur, can be written as: 



huJn — '^^n^n 



sn'^{xi) — sn'^{x2) 



sn{x2)cn{x2)dn{xi)di — sn{xi)cn{xi)dn{x2)d2 K^n. 



One can calculate that hujQ = 0, i.e. \E'o(a;) = ^^0(2^) = e '^^^^Uq, is the lowest 
eigenfunction of H with energy £^0 = 0. Further calculations show that: 



/ 



h{x) 



2(1 - ey, 2V2B{i - 

\u2 ] \ 2V2B; 2 
Diagonalizing the matrix on the r.h.s. and following the algorithm of Subsec- 




tion 4.1, one obtains the eigenfunctions ^n(a5') for (36) 



\l/o = e~'^^^^ dn{xi)dn{x2); 

= e"''^''^ (^cn{xi)cn{x2) + ^ sn{xi)sn{x2)'^ ; 

2 = 6""^''^ (^cn{xi)cn{x2) + ^ sn{xi)sn{x2) 



and their energy eigenvalues: Eq = 0, Ei^2 = 2 — fc^ =f y fc'' + 8-8^(1 — k'^). At 
first sight, this result may be in conflict with the obvious non-negativeness of 
the spectrum of the Hamiltonian H due to its quasi-factorizability (|4]): the 
energy Ei above seems to be negative for some values of B. The resolution 
of this paradox is quite simple, since for the wave functions of interest each 
possible singularity must be normalizable. The wave functions above possess 
the power singularities at {xi ± X2) —>■ via the multiplier e~'^ according 



to expression (37). These singularities are normalizable for B < 1/2, and 
just for only these values of B is the eigenvalue Ei positive. The squared 
wave functions \l/o(a;)2, \l/i(x)2, '$2{xy are given in Fig.lb-ld for the same 
parameter values as in Fig. la. 

The periodicity properties of elliptic Jacobi functions allow us to check that 
all three Bloch eigenfunctions have vanishing quasi-momenta, i.e. they are 
2K-periodic: 

'^nixi+2Kmi, X2+2Km2) = exp {iki{en) ■ 2Kmi + ik2{en) ■ 2X1112)'^ n{xi-, X2), 

where mi,m2 are arbitrary integer numbers, and the quasi-momentum vector 
k{en) = (A;i(e„), A;2(en)) = 0. Although we cannot find the whole structure of 
the spectrum of the model, the zero values of the quasi-momenta hint that 
perhaps the states obtained correspond to the band edges (in analogy with the 
one-dimensional case). In any case, we analytically find three eigenfunctions 



of the generalized two-dimensional Lame Hamiltonian (36) for B^ < 1/2, and 
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one of them, "^oi^x), is certainly the ground state; i.e., the lower edge of the 
bound band. 



4-3 The two-dimensional associated Lame potential 



Since the functional equation ^ is homogeneous, one can consider the more 
general form of (32) and (35), multiplying the functions 0i^2 by an arbitrary 



parameter / and keeping (p^ unchanged: 

01 (x) = 02(a;) = Ik 



2 sn{x)cn{x) 
dn{x) 
cn{\/2x) 
sn{\/2x) 

Thus, the one-dimensional potentials Vi 2 are: 



h+{x) = -0-(x) = B 



(38) 
(39) 



V^i,2(xi,2) = /(/ + l)k\sn\x, 2) + l{l - 1)^2?"^^ - l'k\ 

dn'^[xi^2) 



(40) 



These potentials coincide with the so-called associated Lame potential [T7] . 
with well known band structure [11] , [12] for different (not only integer) values 
of the parameter /. Here we consider onl} the simplest case 1 = 2, where the 
spectrum has two bound bands and the continuous band with the following en- 
ergy values of the band edges and analytical expressions for the corresponding 
wave functions with necessary (anti)periodic conditions under Xi ^ Xi + 2K 
illi: 





= 0; 






= 5- 


3A;2 




= 5- 


2k^ 




= 5- 


2k^ 


64 


= 5- 





Po = Vo = dn^{x)\ 
2v/4^3F; = = ^ {3kHn\x) + a,) 

2V4 - 5A;2 + k^- p, = r^2='^) {?>kHn\x) + A) 

V.= '^,{^kHn\x)+^2) 
^4 = li(3fcW(x) + a2) 



2^4 - + p3 
2^4 - 3F; p4 



1 ' 



where ai,2 = -2 =F a/4 - 3F, /?i,2 = -2 - /c^ =f ^4 - S/c^ + k^. The zero modes 
of q+ are again Un = Pn{xi)T]n{x2). 



This approach seems to be apphcable to the higher values of Z > 2 as well, 
although further calculations will be much more complicated. 
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The two-dimensional potential is: 



'cn\xi) ^ cn^{x2) \ _ ^ 



+ 



dn'^{xi) dii?{x2) 
B{B + V2dn{xi + X2)) B{B - V2dn{xi - X2)) 



(41) 



+ 



2B' 



sn'^{xi + X2) sn'^{xi — X2) 

which is the two-dimensional generalization of the associated Lame potential 
(40). Due to properties of elliptic Jacobi functions, it is again periodic under 
Xj — >■ Xj + 2K ■ rrii with an arbitrary integer mj. 



In order to find a part of spectrum of the system (41), we must study the 
action of the operator h on the functions oOn found above. Again, the matrix 
C is block-diagonal since hujQ = 0, and therefore the lowest energy eigenvalue 
(the lower edge of the first bound band) vanishes. Other zero modes ujn are 
mixed with each other by the action of h: 



h{x) 



UJ2 
UJ3 



I 2Ei 2V2Ban 2V2Bai2 ^ 
2V2Bbii 2E2 2y/2Bbi2 

2y/2Bb2i 2E3 2y/2Bb22 

2^/2Ba2l 2y/2Ba22 2E^ 



where the coefficients ajj, bij have a rather involved form: 



aj(a2 — Oil) 



■1> 



ai 



UJ2 
Us 

3-j . 



(42) 



Li = {l-k')af -2k\3 + ai). 



Diagonalizing the matrix on the r.h.s. of (42), one obtains the energy eigen- 
values En] n = 1,2,3,4 and eigenfunctions "^nix) for (41) analogously to the 
previous Subsection. In this case also the quasi-momenta k{en) for n = 1, 2, 3, 4 
are zero. 



4-4 The two-dimensional trigonometric Razavy potential 



Although the two models studied in the previous Subsections lead to very dif- 
ferent potentials - the two-dimensional Lame and associated Lame potentials- 
the forms of the initial solutions 0i^2 are very similar to each other, and the 
functions (f)± simply coincide. The difference can be described by the param- 
eter 7: 

M-)=-,k^ '"^f^'"^f\ (43) 

dn{x\k) 
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where 7 = 1 for the Lame system, and 7 = 2 for the associated Lame. 
In this Subsection, we consider the hmiting case when 

k^O, (44) 



/3 being a new arbitrary finite parameter, and 0± being the same as in (35), 
but with k 0. Taking into account that in this hmit 

sn{x\k) sinx; cn{x\k) cosx; dn{x\k) — > 1, 

we obtain: 

(^^2(x) = ^sin(2x); (f)+{x) = -(f)^ix) = Bcot{V2x), (45) 
and the one-dimensional potentials Vi 2(3;) takes the form: 

Vi,2{x) = y (1 - cos 4a;) - 2/5 cos 2a;. (46) 



These potentials coincide with the periodic potentials used by M. Razavy 
for the description of torsional oscillations of certain molecules. One must 
choose ^ = —2/5; n = in Eq.(4) of [19] in order to identify the trigonometric 
Razavy potential with our Eg. ([46)). The trigonometric Razavy potential admits 
[21] partial solvability: a few of its band edge levels (eigenvalues and wave 
functions) were found analytically for different values of the integer parameter 
n. Actually, in our case of n = only the lowest band edge eigenfunction with 
eo = was found: 

?7o(a;) = po(a;) = exp(-^cos2a;). (47) 
Our choice for 2 ± leads to the following two-dimensional periodic poten- 



tial (24), which can be considered as a generalization of the one-dimensional 



trigonometric Razavy potential: 
/52 

V{x) = — ^(cos4a;i + cos4a;2) — 2/5(cos2xi + cos2a;2) -l- 

B{B + V2) B{B-V2) (4^) 
sin2(^x+) sm^{V2x^) 

For this potential, the method of SUSY-separation of variables provides the 
analytical expression for the lowest {Eq = 0) band edge wave function: 



^'o(^) =exp {K{x))r]o{xi)po{x2 

B 

smixi — X2' 



sin(a;i + X2) 



v-2 (3 (49) 

exp (— — (cos2xi + cos 2x2)), 
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where k{x) was calculated directly from (26), and rjo, po - from (47). 



5 Discussions and conclusions 



In Section 3 two different symmetries (5*1, 5*2) of the functional equation ([s]) 
were presented. Here we apply them to the Lame and associated Lame poten- 
tials discussed above. 

First, 5*1 symmetry applied to any solution {(pi, 02, 0+, 0-) leads to the poten- 
tial 

which differs from ^ only in the change of variables xi^2 = x±,. Since the 
Laplace operator has the same form in xi^2 coordinates, new wave functions 



'^^^^ of (50) are obtained from the old ones of ^ simply as \I^^'^^^(a;i, X2) = 

In contrast to Si, the effect of S2 symmetry could in principle be more promis- 
ing for the building of new models. Performing the 5'2-transformation of the 



associated Lame potential, generated by (38)-(39) with / = 2, we obtain a new 
potential: 



dn^[xi) dn^[X2) 



-8e + 



^ B{B - V2dn{xi+X2)) ^ B{B + ^/2dn{xi - X2)) _ ^^2 

ClT?{Xi + X2) CTn?{xi — X2) 



(51) 



It is connected to (41) as: 



V^^^\xi, X2; B) = V{xi, X2 + K- -^7|^), 

and the wave functions of V^^'^^ can be obtained easily from those of V by the 
same change of coordinates and parameters. The analogous conclusion is also 
suitable for the Lame potential, the case 1 = 1,. 

In the main text - Section 4 - we in fact considered only one of the superpart- 
ner Hamiltonians H, H from the basic intertwining relations (l|. It is well 
known [Ij that in one-dimensional SUSY Quantum Mechanics these super- 
partners are usually almost isospectral, i.e. they have the same spectra but up 
to normalizable zero modes of the supercharges Q^. The absence of such zero 
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modes in the case of non-periodic potentials means the spontaneous breaking 
of SUSY. The situation is different [TU] for some one- dimensional periodic po- 
tentials, where the two partner potentials V, V are related by a discrete sym- 
metry, but SUSY is not broken. These potentials were called " self-isospectral" , 



and the particular examples are given by the one-dimensional Lame ( 33 ) and 



associated Lame (40) potentials. Another - non-self-isospectral - class of one- 
dimensional periodic models also exists [11], [12], where (quite the contrary) 
the SUSY intertwining relations provide a variety of new solvable periodic 
potentials. 

In the non-periodic models on the whole plane [3], [1], [5] both second-order 
supercharges may have zero modes. It is interesting to consider the two- 
dimensional periodic supersymmetric models from the point of view of 
their self-isospectrality. For example, the superpartner of the generalized Lame 



potential (36) has the form: 

V{x) =2k^{sn'^{xi) + sn^{x2) - 1) + 

^ B{B -V2dn{xi + X2)) ^ B{B + V2dn{xi - X2)) _ 

sn'^{xi + X2) sn'^{xi — X2) ' 

i.e., it differs from V{x) only by the signs in front of the functions dn{xi ±0:2) 
in the nominators. It is easy to check that the reflection {xi,X2) —>■ (xi, —X2) 
merely turns V into V and vice versa: 

V{Xi,X2) = V{xi, -X2). 

Therefore, the spectra of the Hamiltonians H and H coincide, and the two- 



dimensional generalized Lame potential (36) obeys the property of self-iso- 



spectrality. The analogous proof also works for the two-dimensional associated 



Lame potential (41) 



It is significant that while in a certain sense the self-isospectrality for one- 
dimensional models (with first-order supercharges) renders supersymmetry 
useless, this is not the case in the two-dimensional situation (with second or- 
der supercharges). Indeed, since all the two-dimensional models considered 
above satisfy the supersymmetric intertwining relations ([T]) with second-order 
supercharges Q^, the corresponding Hamiltonians H (irrespective of prop- 
erties of the superpartners H) commute with the operators 

R = Q-Q+, 

where are given by This means that all these systems are completely 
integrable - R plays the role of the symmetry operators. It is clear that these 
symmetry operators R annihilate the wave functions constructed in Section 4, 
since they were built simply as linear combinations of zero modes of Q~^. How- 
ever, the action of the operators R on other (as yet unknown) wave functions 
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of H may be very nontrivial. 



In this paper the one- dimensional systems with a finite number of bands were 
represented by the band edge wave functions only. The wave functions in- 
side the bands for the Lame potential, however, are also known. Two linearly 
independent wave functions with energy e are given by: 



H{x ± a) 
Q{x) 



(52) 



where the Jacobi theta-functions, H, B, and the Jacobi zeta-function, Z, are 
defined in the theory of elliptic functions (see [H]). One could use these wave 
functions, instead of band edge functions above, to construct the additional 
eigenvalues of the two-dimensional periodic models. This task seems to be 
much more difficult technically and will be considered elsewhere. Here we 



wish to illustrate how the wave functions (52) coincide with the band edge 
wave functions (34) in the limits e — 0, (1 — k"^), 1. Indeed, these limits 
correspond to the following values of a : K + iK', K, 0, where K' is the 
associated complete elliptic integral K'{k) = K{k') = K{^/l — k?). One then 
has to substitute these limiting values into if, O, Z in (52). As a result, just 



Eqs.(34) are derived. 



Let us also mention the limit k 
ated Lame systems 
1) = sechx 
lead to 



dn{x\k 



1 of the two-dimensional Lame and associ- 
and (41 ). In this limit, sn{x\k = 1) = tanhx; cn{x\k = 
1) = sechx. Therefore, for k 1 Eqs.(38) and (39) 



^'1,2 (a;) = / ■ tanh x; 



where / = 1 for the Lame and / 



(X = - 



(/)_(x) 



B 



sm 

2 for the associated Lame systems. Substitu- 



tion of these expressions into two-dimensional potentials Eqs.(36), (41) gives 



exactly the potential of the two-dimensional generalization of Poschl- Teller 
model. This was studied in detail in [7j, [S], where its partial solvability and 
complete integrability were demonstrated. 



We should notice in conclusion that even in one-dimensional quantum me- 
chanics very limited number of exactly solvable periodic problems is known 
[To] . pT] . [T2] . There is no need to stress the interest of finding analytically 
solvable higher dimensional models like those described in this paper. Such 
models would be very desirable both as a basis for perturbation theory and 
for further study of general properties. The importance of the investigation of 
this kind of model is increasing owing to the development of modern physical 
technologies, which have led to the manufacture of new materials: a variety 
of superlattices, films, quantum two-dimensional dots etc. A study of these 
materials and the corresponding devices should be based on two-dimensional 
(and three-dimensional) Schrodinger equations with periodic potentials. 
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